On the frequency of oscillations in the pair plasma generated by a strong electric 

field 



A. Benedetti, W.-B. Han, R. Ruffini, G.V. Vereshchagin 

ICRANet, P.le della Repubblica 10, 65100 Pescara, Italy, 
ICRA and University of Rome "Sapienza" , Physics Department, 
P.le A. Moro 5, 00185 Rome, Italy. 



o 



Oh 
I 



Abstract 

We study the frequency of the plasma oscihations of electron-positron pairs created by the vacuum polarization in an 
.uniform electric field with strength E in the range 0.2 Ec < E < 10 Ec- Following the approach adopted in [l[ we work 
'out one second order ordinary differential equation for a variable related to the velocity from which we can recover the 
classical plasma oscillation equation when E ^ 0. Thereby, we focus our attention on its evolution in time studying how 
this oscillation frequency approaches the plasma frequency. The time-scale needed to approach to the plasma frequency 
,and the power spectrum of these oscillations are computed. The characteristic frequency of the power spectrum is 
determined uniquely from the initial value of the electric field strength. The effects of plasma degeneracy and pair 
^annihilation are discussed. 
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] The electron-positron pair production in a strong elec- 
tric field is one of the most popular topics in relativis- 
[tic field theory 0]. It begun with the pioneer works by 
Sauter Heisenberg and Euler j^], and by Schwinger 
This effect acquires particular importance when the 
electric field strength E is larger than the critical value 
[Ec = m? / [eh)] such a strong electric field can be reached 
in astrophysical environments, near quark stars 0-0 
neutron stars Strong electric fields up to several 

■percents of the critical value will be reached by advanced 
laser technologies in laboratory jscperiments 
■ray free electron laser facilities [13], optical high-intensity 
laser facilities such as Vulcan or the Extreme Light In- 
frastructure lj|, for a recent review see 15|. Electron 



beam-laser interactions seem also promising in reaching 
high Lorcntz transformed electromagnetic fields capable 
for multiple pair production [l6l |. 

It has been shown that, due to back reaction and screen- 
ing effects of e+e~ pairs on external electric fields, positrons 
'and electrons move back and forth coherently with alter- 
nating electric field: the so called plasma oscillations. In 
[l[ it was pointed out that this phenomenon occurs also 
when E < Ec giving emphasis on the fact that, for over- 
critical (undercritical) field, a large (small) fraction of the 
initial electromagnetic energy is converted into the rest 
mass of pairs, whereas a small (large) fraction is converted 



into kinetic energy. In IT] the case of spatially inhomoge- 
neous electric field has been considered, the emitted radia- 
tion spectrum far from the oscillation region was obtained, 
presenting a narrow feature. 

In this Letter we return to basic equations describing 
pair creation and plasma oscillations in uniform unbound 
electric field. We first derive a master equation for a new 
variable constructed from hydrodynamic velocity, which 
turns out to be second order ordinary differential equa- 
tion. This equation is reduced to the classic plasma oscil- 
lations equation describing Langmuir waves in the limit of 
small electric field. The frequency of oscillations is then 
shown to be almost equal to the plasma frequency, which 
is strongly time dependent in the case under considera- 
tion. Finally, the spectrum of bremsstrahlung radiation 
is computed following 17] and its characteristic feature is 
identified as a function of initial value of the electric field 
strength. 

As in [l| we apply an approach based on continuity, 
energy-momentum conservation and Maxwell equations in 
order to account for the back reaction of the created pairs 
focusing on the range 0.2 Ec < E < 10 Ec. 

We assume that electrons and positrons are created 
at rest in pairs, due to vacuum polarization in uniform 
electric field [l|-[3> [Hl-dQ,] with the average ratc0 per 
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unit volume V and per unit time t 
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where F^^'^ is electromagnetic field tensor, m is electron 
mass. 

This formula is derived for uniform constant in time 
electric field. However, it still can be used for slowly time- 
varying electric field provided the inverse adiabaticity pa- 
rameter fl9'|-[23'l is much larger than one. 



V = fT— = TEpeak > 1, 

UJ Er 



(3) 



where uj is the frequency of oscillations, T — m/uj is dimen- 
sionless period of oscillations. Eq. (jSj implies that time 
variation of the electric field is much slower than the rate 
of pair production. In two limiting cases considered in this 
Letter, Eio = lOEc and i5o.2 — 0-2Ec, we find respectively 
for the first oscillation 7710 ~ 334, and 770.2 = 4.8 x 10^, 
whereas for the last one 7710 — 1.5, and 770.2 — 4. 

Following [1] the conservation laws and Maxwell equa- 
tions written for electrons, positrons and electromagnetic 
field are 



dF'"' 



= -47rJ'', 



(4) 
(5) 
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where n is the comoving number density of electrons, T'^'^ 
is energy-momentum tensor of electrons and positrons 
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where J'' = J^l . 



•^poi total four-current density, 



[/'^ is four velocity respectively of positrons and electrons. 
Electrons and positrons move along the electric field lines 
in opposite directions. 

It has been shown in [ii] that in a uniform electric field, 
from the system ([I])-®, the following system of four cou- 
pled ordinary differential equations may be obtained 



dn 

di 
dp 

di 
dp 

£ 

dE 

di 



fiEv + ^S, 
fiE 4- 77) S*, 

— OTTO! 777; H — 

\ E 



(8) 
(9) 
(10) 

(11) 



where 77 — m^fi is dimensionless number density normal 
ized by the Compton length Ac — l/m, p 



density of positrontQ, p — in^p is momentum density of 
positrons, E = EcE is electric field strength, and t = m^^i 
is time, normalized by the Compton time = l/m. The 

rate of pair production is = ^^iJ^exp {^—-^i velocity 

1^2 

is C = PI p and Lorentz factor is 7 = (l — 7}^) 
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Figure 1: Numerical solution of ([TE| for Eq = 2Ec. This 
figure shows damped oscillations with a frequency increas- 
ing in time. 
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Figure 2: Ratio between number density and maximum 
achievable number density n„iax — l/(47rck) A~'^ for Eq = 
2Ec] it becomes close to unity after 10^ tc- From Fig. [T]and 
Eq. (ITSI) we see that the maxima of u correspond to i? = 
and quenching of pair creation giving rise to fiattening of 77. 
This happens for each oscillation, but it is more evident at 
the beginning due to the double logarithmic scale on this 
figure. 

For the system ([5|)- (fTT|) there exist two integrals (con- 
servation laws) 
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so the particle energy density vanishes for initial value of 
the electric field Eq. Combining together the previous two 



777 "^p is energy 



Total energy density of electrons and positrons is twice this value. 
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equations, we get for the maximum number density of pairs 
that can be created 



(14) 



In [1[ the system (I5|)-([TT|) was reduced to two equations 
using (fT^ -([T ^ and analyzed on the phase plane {E,v). 

Notice that the equation of motion of single particle in 
our approximation is just 



(15) 



where we have defined u — du/dt and introduced a new 
variable constructed from hydrodynamic velocity as ii = 
= "fv/c. Then this equation can be combined with (|lip 
to obtain a single master equation 
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The key point of our treatment is the physical interpre- 
tation of Eq. (1161) which can be rewritten symbolically 
as 

ii + ujlu + ku = 0. (17) 



With constant coefficients Eq. ([T7)) would describe damped 
harmonic oscillations with frequency ujp and friction k. In 
our case uip and k are time dependent, but Eq. (fT6)) still 
possesses an oscillating behavior with damping. With our 
definitions the number density of pairs is 



167ra(l + u2)i/2- 
We then identify ujp in (IT7|) as 



' Sna n 
(l + u2)i/2^ 



(18) 



(19) 



i.e. the relativistic plasma frequencjU. 
The function 



— (1 + ^2) exp 



(20) 



in Eq. ([TT]) accounts for the rate of pair production ([T]). It 
describes the increase of inertia of electron-positron pairs 
due to increase of their number and causes decrease of the 
amplitude of oscillations. 

For small electric fields Eq. dTB]) is reduced to classical 
plasma oscillations equation describing Langmuir waves, 
since in that case k is exponentially suppressed. For this 
reason we expect that as the amplitude of oscillations of 



^The factor Stt is in this formula due to the presence of two charge 
carriers with the same mass - electrons and positrons. This is differ- 
ent from the classical electron-ion plasma where only electron com- 
ponent oscillates and the corresponding factor is twice smaller. 



electric field gets smaller the frequency of oscillations uj 
tends to the plasma frequency Wp. 

In [l[ we solved numerically the system of four coupled 
ordinary differential equations (I5|- (|lip . Now it is possible 
to solve just one second order differential equation ([TCI) 
which allows us to study its asymptotic behavior as well. 

We solve numerically Eq. (flB)) with the initial condi- 
tions {t(0) = Eq and u{0) — 0, corresponding to no pairs in 
the initial moment, taking for initial electric field strength 
Eq = {0.2,0.5,1,2,5,10}. Once this equation has been 
solved, we have the solution for the number density from 
(fT8|l and for the plasma frequency by means of (fT9)) . In 
Fig. [l]we show the evolution in time of u where the am- 
plitude of the oscillations decreases, while its frequency 
increases with time. 

We present the number density of pairs in Fig. [2]for the 
case _Bo = 2 i^c as a fraction of the maximum achievable 
value rimax- In all the cases under interest, this number 
is asymptotically achieved indicating that the final result 
of the process will be the complete conversion of the elec- 
tromagnetic energy density into the rest mass of the pairs. 
Moreover, looking at Fig. [2]we recover the result obtained 
in [l|; in fact we can see that after the first oscillation, 
higher is Eq larger is n. This means that the first oscilla- 
tion gives the leading contribution to the process in which 
the electromagnetic energy of the field is converted in the 
rest mass of pairs, with a moderate contribution to their 
kinetic energy for Eq > Ec- The values of the half pe- 
riods of the first oscillation for each considered case are 
represented in Fig. El by triangles. 

We computed the frequency of the i-th oscillation as 
tj* — 1^1/2^ where is the corresponding half period, 
calculated considering the time interval between the i-th 
and (i -I- l)-th subsequent roots of u, see Fig. [1] 

Notice that ujp is an oscillating function of time, due to 
the presence of u in (fTO)) : besides the frequency of these os- 
cillations increases in time. For this reason, in order to get 
a smooth function we calculated the average of Cjp . We use 
this new function a)"" for the plasma frequency to make a 
comparison with the frequency of oscillations of pairs. In 
Fig. 131 for the case Eq — 2Ec, the blue area represents the 
plasma frequency as defined by (jl9p . the yellow curve is 
its average w™, while the pairs oscillation frequency w is 
represented by the red curve. For the same initial electric 
field, in Fig. [S]the trend of the ratio w/cjp" is shown, which 
indicates that the averaged plasma frequency is achieved 
asymptotically as expected from (ITBl) . Notice that the os- 
cillation frequency a) is always smaller than a)"" since the 
number density of pairs is constantly increasing with time 
during each oscillation cycle. 

We computed the power spectrum of radiation in the 
far zone, assuming dipole radiation following jl7| . The 



power spectrum, namely the energy radiated per unit solid 
angle per frequency interval and per unit volume, is given 

by 

P(^)^^^2a|D(^)P, (21) 
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where the ampUtude D(a}) is proportional to the Fourier 



transform of the electric current time derivative [17 1 
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The electric current is simply related to the new variables 
by the following expression 



J = 2 ^/ah ■ 



(1 + ^2) 



r,2\i/2- 



(23) 



From the Fig. [5] it is clear that the main contribution is 
given by the final and fastest oscillations which last for a 
longer time. Therefore, the power spectrum shows a peak 
close to the plasma frequency being reached asymptoti- 
cally. We can easily estimate the frequency corresponding 
to this peak combining Eq. dTi|) with Eq. as 



^peak 



Eq 

71 



(24) 



with the corresponding energy hujpeak — 0.72 Eq/Ec MeV. 

The energy loss due to the dipole radiation for E — 2Ec 
considered in Fig. [5] for t = lO^tc is less than one percent. 

Once we know the frequency of the pairs oscillations 
and the plasma frequency, we obtain their ratio as it is 
shown in Fig. [5j Besides, we use w/Wp" in order to com- 
pute the characteristic time scale ta needed for the pairs 
oscillation frequency to reach the plasma frequency. This 
has been done considering the ratio between w/a;™ and 
its time derivative; the result of this procedure gives us a 
numerical function from which we have taken the average. 
For all the considered cases, this quantity is shown in Fig. 
[31 from which we understand that the general trend is that 
larger is the initial electric field, larger will be the starting 
oscillation frequency. 

It is worth noting the effect of degeneracy on the pair 
production. One may think that when concentration of 
pairs reaches the maximum allowed value by the Pauli 
principle the pair production is blocked. For particles at 
rest this would happen when two pairs with opposite spins 
occupy a Conipton volume. Considering asymptotic num- 
ber of pairs given by (jl4p one finds that it would happen for 
E > Ay/oTrEc — O.GEc- However, one has to keep in mind 
that particles produced at rest are accelerated by external 
electric field and thus leave the quantum state with zero 
momentum which can be subsequently filled by a new pair. 
These effects are independent since they operate in ortog- 
onal directions of the phase space, so one can estimate the 
value of external electric field at which phase space block- 
ing occurs by comparing their rates. Such analysis gives 
us the following inequality 



1 ~ / TT 

—^E exp — ^ 
47r3 \ E 



> 1, 



(25) 




Figure 3: In this figure are shown, for all the cases under 
interest, the half period of the first oscillation ti (trian- 
gles), the characteristic time-scales ta (squares) needed to 
the pairs oscillation frequency to reach the plasma fre- 
quency and the time t^ (circles) which satisfies the condi- 
tion T{tj) ~ 1 for the optical depth defined in (1^^ . The 
value of ta for the case Eq — 2Ec is shown in fig. [5] by the 
vertical line. 
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Figure 4: In this plot the blue area represents the plasma 
frequency as defined in Eq. (|19p: it appears like a contin- 
uum because of the fast oscillations. The yellow curve is 
its average in time w"" which can be compared with the 
pairs oscillation frequency w given by the red curve. This 
plot corresponds to the case Eq — 2Ec- 
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Figure 5: Behavior of the ratio Cj/Cj^" in tinie for Eq = 
2Ec- The plasma frequency is attained asymptotically be- 
cause of the limit Sit — ^ oo) — >■ 0. The vertical line cor- 
responds to the time-scale needed to attain the plasma 
frequency in this specific case. 



having the solution E > 127Ec, which is much higher than 
electric fields considered in this Letter. 

Another effect, relevant for large enough electric field, 
is interaction of pairs with photons discussed in some de- 
tails in [l|, HSl- '^^^ estimate the optical depth for 
electron-positron annihilation as 



r{t) 



(Tt 

7 



nvdt 



STra^ 



(l + u2)3/2 



fi dt. 



(26) 



where ctt is the Thomson's cross section, and we approx- 
imated a ~ Equating (|26p to unity we find the 
timescale tj at which the probability of electron to interact 
with positron and create a pair of photons reaches unity. 
From that time moment interaction of pairs with photons 
can no longer be neglected. This timescale is represented 
in Fig. [3] by circles. 

Summarizing the informations presented in Fig. [3] we 
conclude that independent on the initial value of the elec- 
tric field, there is a hierarchy between the following time 
scales ti < t-y < ta- It means that many oscillations occur 
before electron-positron collisions turn out to be impor- 
tant, thus justifying our coUisionless approximation. The 
condition t-y < ta means that the estimation of the max- 
imal frequency of oscillations is an approximate one: 
photons produced by interaction of pairs will also distort 
the spectrum shown in Fig. [51 

As long as electric field does not reach critical val- 
ues for creation of muons and pions their production from 
electron-positron collisions 25]-[26] is suppressed because 
of two different mechanisms. Both these processes have 
a kinematic threshold given by the rest mass of the pro- 
duced particles. For this reason the Lorentz factor of the 
relative motion of colliding electron and positron should 
exceed ^ 10^, restricting initial electric fields to be under- 



hand, the number density of pairs produced is exponen- 
tially suppressed for undercritical fields. Besides the cross 
section for all these processes decreases as cr oc 7"^ which 
further decreases the rate of electron-positron collisions. 

It has been claimed recently [l^l that critical Schwinger 
field could never be reached in high power lasers due to oc- 
curence of avalanche-like QED cascade operating mainly 
through nonlinear Compton scattering combined with non- 
linear Breit- Wheeler process 28l 
the trident process 



see 



also 0, and via 
As soon as one single pair is 
generated by the Schwinger process such electromagnetic 
cascade of secondary electron-positron pairs is expected to 
deplete the electromagnetic energy thus preventing further 
pair production from vacuum. The requirements for the 
avalanche to occur are twofold: a) the probability to emit 
photon should not be suppressed and b) the photon must 
be energetic enough to produce pair by interaction with 
another photon. It is shown that for a s pec ific electromag- 
netic field configuration considered in (28l [29| as well as 
namely circularly polarized standing electromag- 



27| 



netic wave both these conditions may fulfill for undercriti- 
cal electric field E < Ec- However, as it was shown in (3^ 
for linearly polarized standing wave such electromagnetic 
cascade is not expected to dominate over the Schwinger 
process. It is easy to understand these results looking at 
the energy loss rate of charged particle in classical electro- 
dynamics 
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3 TO 
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(E + V X H)^ - (E- v) 



(27) 



When magnetic field is absent (the case considered in [28|, 
[29! and 27]) if directions of particle velocities and electric 
field are collinear the radiation loss turns out independent 
on particle energy. In such case, as we have shown pre- 
viously [1] for overcritical electric field the radiation loss 
is smaller than the rate of energy conversion from elec- 
tromagnetic field to electron-positron pairs via Schwinger 
process. Notice that in the case of plasma oscillations con- 
sidered in this Letter the velocity and acceleration vectors 
are indeed collinear, so curvature radiation considered in 
^28. . J9] does not occur. When electric field changes with 
time not only its amplitude but also direction, as for in- 
stance in circularly polarized electromagnetic wave, the 
acceleration and velocity vectors become misaligned and 
curvature radiation becomes much more efficient due to 
quadratic dependence on particle energy in ([?f)) . We also 
notice that the backreaction of electron-positron pairs on 
the initial electric field, which is the top ic of the present 



and [27| 



see 



Letter, is not taken into account in 
however [sst . 

To summarize, the study of plasma oscillations due to 
the vacuum polarization in uniform electric field can be 
reduced to the analysis of a single second order ordinary 
differential equation for the variable constructed from hy- 
drodynamic velocity ii = "fv/c. All the other physical 



critical, Eq < Ec, see Fig. 3 in Ref. On the other quantities of interest can be obtained from the solution of 
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Figure 6: The power spectrum in arbitrary units has been 
obtained using Eq. (HI]). The peak almost corresponds to 
the maximum plasma frequency (jl9l) corresponding to the 
maximum achievable number density (|14p. indicated here 
by the vertical line. The main contribution is given by the 
oscillations which last for a long time, namely when the 
asymptotic limit of the plasma frequency is attained. 



(|16|) . This reduction allows to study the evolution of the 
system for a long time. 

As expected, the plasma frequency is reached asymp- 
totically for all the considered cases 0.2 Ec < Eg < 10 Ec- 
The difference between them rely on the time scale the sys- 
tem needs to approach the plasma frequency as it is shown 
in Fig. [31 In particular, larger is the initial electric field, 
shorter will be the time scale to get the plasma frequency. 

Surprisingly we find that, for all the cases we have con- 
sidered, Cj ~ uip" even for the very first oscillations, when 
we are far from the asymptotic case S{i oo) — > when 
we expect w = w^" from the analysis of Eq. ((TB| . 

The characteristic feature of the power spectrum of 
dipole radiation occuring due to plasma oscillations is shown 
to be located close, but always below, to the plasma fre- 
quency. The left tail in Fig. |6]is due to the first oscillations 
with frequencies smaller than Wp", while the main contri- 
bution is due to the final evolution when the pairs oscillate 
almost with the same frequency close to Wp". 

The upper limit to the optical depth to pair annihila- 
tion into photons is obtained, showing that it never exceeds 
unity for E < AbE^. 
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